The thermodynamical liquid-glass transition in a Lennard-Jones binary mixture 
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We use the results derived in the framework of the rephca approach to study the liquid-glass 
thermodynamic transition. The main results are derived without using replicas and applied to the 
study of the Lennard-Jones binary mixture introduced by Kob and Andersen. We find that there 
is a phase transition due to the entropy crisis. We compute both analytically and numerically the 
value of the phase transition point Tk and the specific heat in the low temperature phase. 



In recent times there have been many progresses in the 
analytic understanding of the thermodynamics transition 
of glasses [l]-0 . Beyond technicalities the basic assump- 
tion is that a glass is very near to a frozen liquid. This 
quite old statement can be rephrased as a liquid is very 
near to a heated glass. In other words the configurations 
of a glass at low temperature are not far from those of 
a liquid. Therefore, if we use some smart method to ex- 
plore the phase space in the liquid phase, we can find the 
properties of the low temperature glassy phase. 

This strategy has been put into action using the replica 
theory. One finds, as output of an explicit computation, 
that the glass transition is characterized by the vanishing 
of the configurational entropy (the so called complexity) 
and the low temperature phase is described by one step 
replica symmetry breaking |5[| with a non- vanishing non- 
ergodic parameter at the phase transition point. 

Replica theory is a very powerful tool, but it has the 
disadvantage that many of the underlying physical hy- 
pothesis cannot be seen in a clear way. We will rederive 
some of the main result of ||l| without using the replica 
formalism [pj. We suppose that below some tempera- 
ture (to be identified with the mode coupling transition 
M ) , the phase space of the system can be approximately 
divided into regions (which we will call valleys) which 
are separated by high barriers [raM| . This approximation 
becomes better and better when we decrease the tem- 
perature and it becomes exact below the thermodynam- 
ics phase transition temperature Tk of the glass [10|, 
where the viscosity should diverges pd| (unfortunately 
no quantitative predictions have been obtained micro- 
scopically on the behaviour of the viscosity beyond the 
Adams Gibbs argument). In a first approximation each 
valley can be associated to one inherent structure, i.e. 
one minimum of the potential energy ||8|,|l2|,^ . 

Let us consider a generic system with N particles with 
Hamiltonian H{C), C denoting the generic configuration 
of the system. The partition function both in the liquid 
phase at low temperature and in the glass phase can be 
written as 



Z(/3)=^exp(-/37V/„ (/?)) = 

a 

dM{f,P)eM~PNf), 



(1) 



where fa (/3) is the free energy density of the valley labeled 
by a at the temperature (3~^ and N{f, (3) is the number of 
valleys with free energy density (per particle) less than /. 
In the glassy phase this sum is dominated by the valleys 
with minimal free energy, while in the liquid phase an 
exponentially large number of valleys contributes to the 
partition function. It is normally assumed that in glass 
forming systems, for large N and / > foiP), we have 
that 



AA(/,/3)=exp(7VI](/,/3)), 



(2) 



where the configurational entropy, or complexity, S(/, (3) 
is positive in this region and vanishes at f — fo{f3). The 
quantity fa{(3) is the minimum value of the free energy: 
AA(/,/3)iszerofor/</o(/3). 

The key point is the computation of the function 
!](/, /3) for / > /o(/3). It can be done using liquid the- 
ory because in this region we are still in the liquid phase. 
The location of the zero of S(/, /3) will tell us the value 
of the free energy in the glassy phase. 

At this end it is convenient to consider the generalized 
partition function 

Z{r, 13) = exp(-iV7$(7; /3)) = ^ exp(-7iV/, (/?)). (3) 

a 

The physical meaning of 7 is will be clearer later. 
It is evident that 

,^r,P)-,f~m.a / = ^M^2iM. (4) 

The complexity can be simply obtained from <I>(7;/3) in 
the same way as the entropy can be obtained from the 
usual free energy. 

The crucial step consists in writing 



J dCexp {-jH{C) - NjfiP, C) + iV7/(7, C)) = 

dC exp (-iH{C) - iV7/(/3, C) + ^¥7/(7, C)) , (5) 

where /(^, C) = /(/3, C)-/(oo, C) and /(/3, C) is a func- 
tion that is constant in each valley and it is equal to the 
free energy density of the valley to which the configura- 
tion C belongs. In deriving eq. (ph we have assumed that 
the value of 7 is high enough that all the configurations 
C, which contribute to the integral for large N , belong 
to some valley: We have also assumed that there is a 
one to one correspondence among the valley at inverse 
temperatures (3 and at 7. 

Before entering into the computation of /(/3, C) it is 
useful to make the so called quenched approximation, i.e. 
to make the following approximation inside the previous 
integral: 



exp(-^/(/3,C))=exp(-A^(/3)), 



(6) 



where f-y{P) is the expectation value of f{P,C) 
taken with the probability distribution proportional to 
exp(— 7ff(C)). The quenched approximation would be 
exact if the temperature dependance of the energy of 
all the valleys would be the same, apart from an over- 
all shift at zero temperature. In other words we assume 
that the minima of the free energy have different val- 
ues of the free energy but similar shapes. More refined 
computations can be done to compute systematically the 
corrections to the quenched approximation. This approx- 
imation would be certainly bad if we were using the free 
energy /(/3,C) at the place of /(/3, C) because the zero 
temperature energy strongly varies when we change the 



minimum. 



We finally find 



$(7;/3)=Fi(7) + /^(/3)^^(7), 



(7) 



Fl{i) being the free energy of the liquid {Sl{i) will be 
the entropy of the liquid) . A simple algebra shows that 



E(7; (3) = Sl{i) ~ S,{j) + f'{j) - fjp), 



(8) 



where/;(/3)=a^(/3)/97- 

In the liquid phase, we find out that the configurational 
entropy is given by 



E(/3) = S(/3;/3)=5l(/3)-5M/3), 



(9) 



as expected: the entropy of the liquid is the entropy of 
the typical valley plus the configurational entropy. 

The thermodynamic transition is characterized by the 
condition 



In the glassy phase the free energy can be found by 
first computing the value of 7(/3) such that 



S](7(/3);/3)-0 



fll) 



S(/3. 



K 



0. 



(10) 



and evaluating the corresponding free energy. 

The quantity 7(/3) is the inverse of the effective tem- 
perature of the valley. It is easy to show (following ||] ) 
that the previous formulae are completely equivalent to 
the replica approach. 

A strong simplification happens if we assume that the 
entropy of the valley can be evaluated in the harmonic 
approximation where we only keep the vibrational con- 
tributions. In this case we obtain for a system with M 
degrees of freedom that the harmonic entropy of the val- 
ley near to a configuration C is given by 

S{P{C)) = y In {^^ ^Tr (ln(H(C))) , (12) 

where 7i(C) is an M x M Hessian matrix (e.g. if H 
depends on the coordinates Xi we have that Tii^k ~ 
d'^H/dxidxk)- 

If our approximation were fully consistent, we should 
find that all the eigenvalues of Ti (the so called INN, 
Instantaneous Normal Modes jlj]) are positive. This is 
not the case, however the number of negative eigenval- 
ues becomes very small at low temperature, still in the 
liquid phase, signaling that valleys can be approximately 
defined in this region. 

How can we evaluate the harmonic entropy in this 
case? Two rather similar possibilities are: (a) We com- 
pute Tr(|ln(7^(C))|) instead of Tr (ln(7^(C))). (b) We 
find the minimum of the Hamiltonian which is the nearest 
to C (i.e. the corresponding inherent structure) and we 
use the spectrum of the Hessian at this point. We have 
checked numerically that the two methods give rather 
similar results y and we will follow the second one in 
the numerical simulations presented in the paper. 

The framework has been set. We now segue into an 
explicit computation. Two possibilities are open: (a) We 
do an analytic computation of all the quantities which 
appear in the previous equations; (b) We extract them 
from numerical simulations. The first possibility is open 
only in relatively simple systems and further approxima- 
tions are needed, the second one is viable for all systems. 

In these letter we shall explore both possibilities in the 
case of a realistic model for glasses, the binary mixture 
of particles (80% large particles, 20 % smaller particles) 
interacting via a Lennard-Jones potential, introduced by 
Kob and Andersen pq] . This Hamiltonian should mimick 
the behaviour of some metallic glasses and it is one of the 
best studied and simplest Hamiltonian which do not lead 
to crystalization at low temperature. 

The numerical procedure does not present any seri- 
ous difficulty. We have studied via Monte Carlo simu- 
lations a system of TV = 260 particles, in a cubic box 



with periodic boundary conditions at density p = 1.2, 
starting from /3^/* = 0.02 and increasing f3 by steps of 
A/?!/" = 0.02; we performed up to 4 10*^ MC steps at 
each (3. We need to simulate the system at high tem- 
peratures in order to compute the entropy of the hq- 
uid (we use as reference entropy that of a perfect gas 
at /3 = 0). The entropy is obtained using the formula 
S{(3) = 5(0) + ^1^ d/3'{E{P) - E{(3')). Given an equilib- 
rium configuration the near minimum of the potential is 
found by steepest descendent and the computation of the 
780 eigenvalues of 7i(C) does not present any particular 
difhculty. 

The results for the total entropy of the liquid and for 
the harmonic part are shown in fig. (1) as function of 
T""* (a more detailed description of the simulations can 
be found in ref . [Q ) . The entropy of the liquid is remark- 
able linear when plotted versus T~'^, as it happens in 
many cases |l6| . 

In fig (2) we show the configurational entropy as func- 
tion of T~'*. We fit it with a polynomial of second de- 
gree in r^**. The extrapolated configurational entropy 
becomes zero at a temperature Tk = -31 ± .04, where the 
error contains systematic effects due to the extrapolation 
(similar conclusions have been reached in ref. |lj]). 
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FIG. 1. Analytical entropy of the liquid (upper line) com- 
pared with the numerical one, and analytical harmonic en- 
tropy (lower line) compared with the numerical results. The 
horizontal axis is T"*. 

There are many methods to compute analytically the 
free energy in the liquid phase which lead to integral 
equations for the correlation functions. Here we follow 
a simple procedure proposed in ref. [|l8| . It mixes the 
HNC (hypernetted chain) and MSA (mean spherical ap- 
proximation) closures by means of a single parameter 
a{T) that is chosen in order to reduce thermodynamic 
inconsistencies, minimizing the difference between two 
different ways of computing the compressibility. This 
technique allow us to compute the internal energy in the 



liquid with a reasonable approximation. The resulting in- 
tegral equations for the correlations function of the two 
kinds of particles are transformed in a set on non linear 
differential equations by discretizing them with a space 
resolution of 2~^ and a large distance cutoff of 2'^ (we 
have varied these parameters and checked that this choice 
gives a reasonable accuracy). The resulting equations in 
3-512 unknown have been solved using a package from 
the IMSL hbrary. 

The computation of the spectrum is much more in- 
volved. We follow a simple approach which becomes ex- 
act if we assume that the range of the Hessian is much 
larger than the typical interatomic distance and we use 
the superposition approximation 

5(2:1X2X3) = g{xi - X2)g{x2 - X3)g{x3 - xi) (13) 

in the cases where it is needed |H,|10. In order to com- 
pute the eigenvalues of the Hessian we compute analyti- 
cally the moments of the eigenvalue distribution 



Mk = Tr(H'=). 



(14) 



Each moment can be written as the appropriate integral 
over the correlations functions. In this long range ap- 
proximation some terms are dominant over the others. 
We keep only those terms and we introduce the super- 
position approximation into the appropriate places. We 
use the large range expansion to select the diagrams p9| . 
After some computations (which we do not report here 
iljilll]) we find a spectral density p(e) that has support 
only in the region of positive eigenvalues e (real frequen- 
cies to = y/e) and goes to zero as w as expected. 




FIG. 2. Analytical value of the complexity (upper line) 
compared with the numerical one (-1- points), as functions of 

We can now put everything together: the final ana- 
lytic predictions for the liquid and harmonic entropy are 
shown in fig. (1). The liquid entropy turns out to be very 
good, while there is a minor discrepancy for the harmonic 



entropy, likely due to the rather strong approximations 
we have done in the analytic computation. The analytic 
configurational entropy is shown in fig. (2). It becomes 
zero at Tk = .32, which is our analytic prediction for 
the thermodynamic transition (as a check we have fit- 
ted the analytic configurational entropy using the same 
procedure as for the numerical one and we have found 
Tk = .34). 

We have in our hands all the tools to compute ana- 
lytically the free energy in the low temperature case. In 
fig. (3) we show the specific heat coming from the x 
dependent part of the Hamiltonian as a function of the 
temperature (we must add the momentum contribution 
3/2 in order to get the total specific heat). Very similar 
results are obtained by using the extrapolated numerical 
entropies at the place of the analytic ones. We notice 
that the Dulong Petit law is extremely well satisfied in 
the low temperature region. Indeed in the harmonic ap- 
proximation the Dulong Petit law would be exact if we 
neglect the 7 dependence of Sj{(3). The value of ^{(3) 
weakly depends on /?: its value in the limit /3 — > 00 is 
only about 10% higher that its value (i.e. Pk) at the 
transition temperature. 

The results presented here are an explicit realization 
of a transition driven by an entropy crisis, which has 
been firstly implemented microscopically in the random 
energy model ||2^. The specific heat has a jump down- 
ward, when we decrease the temperature, which is the 
opposite of the typical behaviour in transitions charac- 
terized by the onset of a conventional order parameter (in 
mean field theory ferromagnets, superconductors, have a 
jump upward when we decrease the temperature). 
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FIG. 3. The specific heath coming from the x dependent 
part of the Hamiltonian as a function of the temperature 

We stress that it is possible to remove the approxima- 
tion of using the harmonic entropy. There are no serious 
difficulties in computing numerically the true entropy of 
the valleys; this has been done for a binary mixture of 
soft sphere and the results are very near to the harmonic 



ones 1^. Analytically we can also expand in the anhar- 
monicity parameters. It is also possible to take care of 
the fiuctuations of the entropy from valley to valley and 
go beyond the quenched approximation. It is quite likely 
that these effects will not strongly change the results. 
The most important step would be to reach a better the- 
oretical control on the spectrum of the INN in the liquid 
phase, maybe combining the approximations used here 
with the low density expansion of pl| . 

Summarizing we have found a method which is able 
to use liquid theory method in the glasses phase putting 
in practice the old adage a glass is a frozen liquid. We 
are able to compute with a reasonable approximation the 
thermodynamics and with a little more effort we can com- 
pute the static and the dynamic structure functions. 

We owe a lot to M.Mezard and we are very happy 
to thanks him. We are also very grateful to W.Kob 
and F.Sciortino for their suggestions and to A.Cavagna 
and I.Giardina for useful discussions. B.C. would like to 
thank the Physics Department of Rome University 'La 
Sapienza' where this work was partially developed dur- 
ing her PhD. 
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